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We  construct  a  nonhnear  kinetic  equation  and  prove  that  it  is  well-adapted  to  describ 
general  multidimensional  scalar  conservation  laws.  In  particular  we  prove  that  it  $s  well-posei 
uniformly  in  e  -  the  microscopic  scale.  We  also  show  that  the  proposed  kinetic\ equation  is 
equipped  with  a  family  of  kinetic  entropy  functions  -  analogous  to  Boltzmann’s  microscopic 
77-function,  such  that  they  recover  Krushkov-type  entropy  inequality  on  the  macroscopic 
scale.  Finally,  we  prove  by  both  -  BV  compactness  arguments  in  the  multidimensional  case 
and  by  compensated  compactness  arguments  in  the  one-dimensional  case,  that  ihe 


it  ihe 

density  of  kinetic  particles  admits  a  ‘continuum’  limit,  as  it  converges  strongly  withV  jt  0  to 
the  unique  entropy  solution  of  the  corresponding  conservation  law.  /  ■'  j  . 
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1.  INTRODUCTION 


Consider  the  scalar  multi-dimensional  conservation  law 

(L1)  +  J2  ^[^(ufot))]  =  0,  ( x,t)eRi  x  Rf,  Ai(-)eC\ 

with  given  initial  conditions  u(x,  t  =  0)  =  u0(x).  We  are  concerned  here  with  a  Boltzmann- 
like  kinetic  equation  which  describes  (1.1),  as  its  microscopic  scale,  e  >  0,  tends  to  zero. 
To  this  end  we  introduce  a  scalar  function,  /e(x,u,  t),  which  can  be  viewed  as  a  microscopic 
description  for  the  density  of  particles  located  at  (x,t)eR*  x  Rf  with  speed  veR.  Starting 
with  given  initial  distribution,  fe(x,v,  0),  our  kinetic  model  evolves  according  to 

C1-2®)  [dt  +  a(v)  ■  dx]fe(x,v,t)  =  -  fe(x,v,t)\. 

Equation  (1.2a)  tells  us  that  the  particles  are  transported  along 

a(u)  dx  =  J2  a'(v)jrr>  °.(')  =  4(')» 

i=  l 

and  that  their  collisions  are  governed  by  the  nonlinear  kernel  on  the  right.  Here, 

(1.26)  ue(x,t)  =  f  fe(x,  v,  t)dv, 

Jv 

denotes  the  local  density  of  particles  at  a  given  (x,  t )  location,  and  the  ‘equilibrium  function’, 
Xu.(x,t)(u).  is  the  signature  of  ue(x,t),  i.e., 

( i  2c)  Y  (v\  -  {  s8nu>  if  (u  -  v>  >  0. 

K  ‘  Xu{  J  ~\  0,  if  (« -  v)v  <  0  • 

The  classical  example  of  a  kinetic  model  is  of  course  the  Boltzmann  equation  [1].  Equation 
(1.2)  is  closely  related  to  the  B.G.K.  model  of  Boltzmann  equation.  Existence  theory  for 
Boltzmann  equation  and  its  simplified  B.G.K.  model  can  be  found  in  [6],  [11],  respectively. 
In  both  cases,  however,  the  question  of  convergence  of  the  macroscopic  moments  to  weak 
solutions  of  compressible  Euler  equations  is  still  an  open  problem.  (Consult  [3]  regarding 
an  affirmative  answer  to  this  convergence  question  in  the  case  of  strong  solutions.)  In  this 
paper  we  restrict  our  attention  to  the  simpler  scalar  case,  and  we  show  that  the  proposed 
kinetic  equation  (1.2)  is  well  adapted  to  describe  strong  as  well  as  weak  solutions  of  (1.1)  as 
£  1  0. 

The  paper  is  organized  as  follows.  In  Section  2  we  show  that  the  kinetic  equation  (1.2)  is 
well-posed  in  L°°(Rf  \  Ll(R*  x  Rv )).  Next,  we  borrow  our  terminology  from  the  framework 
of  Boltzmann’s  kinetic  equation.  The  microscopic  scale,  £,  in  (1.2)  can  be  viewed  as  the 
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mean  free  path.  In  Section  3  we  prove  that  the  continuum  or  ‘fluid’  limit  of  the  local  density 
of  particles,  lim^o  ue(x,  t)  is  the  unique  entropy  solution  of  (1.1).  A  kinetic  construction  of 
conservative  solutions  was  carried  out  by  Giga  and  Miyakawa  [7].  In  fact  their  construction 
is  nothing  but  a  fractional  splitting  solution  of  our  kinetic  equation  (1.2),  namely,  a  kinetic 
approximation  is  constructed  by  a  succession  of  small  time  steps,  in  which  we  first  transport 
and  then  project  the  particles  distribution  according  to  (1.2).  Here  we  improve  on  [7]  by 
identifying  the  underlying  kinetic  equation  which  corresponds  to  (1.1).  It  is  also  shown  here 
that  this  kinetic  equation  is  equipped  with  (a  family  of)  kinetic  entropy  functions  which  play 
an  analogous  role  to  Boltzmann’s  77-function.  In  particular,  Krushkov  entropy  inequality 
[8],  [9]  is  recovered  in  the  ‘fluid’  limit  e  j.  0. 

In  Section  4  we  revisit  the  question  of  the  ‘fluid’  limit  in  the  case  of  one-dimensional 
kinetic  model.  Here  we  show  that  the  compensated  compactness  theory  of  Tartar  [13]  can 
be  adapted  as  an  alternative  approach  for  providing  an  affirmative  answer  to  the  question 
of  macroscopic  convergence.  The  compensated  compactness  arguments  allow  us  to  pass 
to  the  continuum  limit  with  minimal  L 1  f)L°°  information  about  the  distribution  function, 
fe,  which  may  still  oscillate  around  the  ‘equilibrium  function’  Xu-  Finally,  in  Section  5,  we 
indicate  the  extension  of  our  results  to  the  inhomogeneous  case,  in  the  presence  of  a  (possibly 
stiff)  source  term. 

2.  THE  KINETIC  EQUATION  IS  WELL-POSED 

Let  us  rewrite  (1.2)  in  the  form 

d  ,  ..3.1,1 

(21)  Jj.fc  +  a(u)  •  ^/e  +  -/.  -  "Xu., 

separating  between  its  linear  part  on  the  left  and  its  nonlinear  kernel  on  the  right.  By 
Duhammel’s  principle,  (2.1)  admits  the  following  equivalent  integral  representation 

(2.2)  fe(x,v,t)  =  e"*/e( x  -  ta(v),  v,  0)  +  ^  j  q  e(T_t)/eXu.(x-(t-r)o(v).r)(v)dr. 

The  question  of  existence  of  a  kinetic  solution  of  (1.2)  is  now  transformed  into  that  of  a  fixed 
point  solution  for  the  RHS  of  (2.2).  Fixing  T,T  >  0,  we  seek  such  a  fixed  point  solution 
in  Z,oo([0,7’];  Ll{R*  x  /?„)).  To  this  end,  we  let  /e  and  ge  be  two  different  solutions  of 
(1.2a)  with  corresponding  densities  ue{x,i)  =  fv  fe(x,v,t)dv  and  wt{x,t)  —  Jv  ge(x,v,t)dv. 
By  (2.2),  their  difference  does  not  exceed 

«/.(*,  0  -  9'(x>  v>  OIUmkSxH.)  <  e_t/‘!| /«(*,  V,  0)  -  ge(x,  v,  0)||il(^x^)+ 

+  7  /  e(T_t)/e||x«.(«,r)(v)  -  X«M«,T)(v)|Ui(Ki*Kv)*-- 

c  Jt~ 0 
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Using  the  properties  of  the  signature  function,  x>  we  therefore  conclude 


II  fe(x,v,t)  -  ge(x,u,<)||L1(HixRi()  <  e  </e||/e(x,u,0)  -  ge(x, V, 0)|| Lx (R* x Rv)+ 

(2-3)  f/ 

+  (!  ~  e  )j<r<tll  f‘(X>V’T)  ~  9c(x,V,T)\\Ll(RURv). 

The  inequality  (2.3)  shows  that  the  fixed  point  iterations 

(2.4)  f?^(x,v,t)  =  -  ta(v),v,  0)  +  -  f  e^leX^{x-{t-r)a(v),r)(v)dr, 

E  Jr= 0 

are  contracted  (with  a  contraction  factor  of  1  —  e-t/e)  to  a  fixed-solution  solution  of  (2.2). 
Moreover,  by  (2.3)  this  kinetic  solution  is  unique  and  continuously  dependent  on  the  initial- 
data,  for 

(2-5)  max  ||/e(x,u,t)  -  5e(rc,u,t)||Li(^xflv)  <  ||/e(x,  u,  0)  -  ge(x,  v,  0)||L.(^XiRv). 

We  summarize  this  by  stating 

THEOREM  2.1.  The  kinetic  model  (1.2)  is  well-posed  in  L°°(R( ;  Lx(/2f  x  i?„)).  Moreover, 
its  solution  operator  is  nonexpansive  in  this  topology,  i.e.,  (2.5)  holds. 

We  close  this  section  with  several  remarks. 

1-  ^°° -bound.  To  see  that  the  solution  operator  associated  with  the  kinetic  model  (1.2)  is 
uniformly  bounded,  we  use  (2.2),  obtaining 

ll/e(’>  0IU°°(Hg)  ^  e  </e||/e(')'y>0)IU“(K2)  + 

(2.6)  +  (1  -  e“t/e)  max  ||xu.(..r)(v)||x,«>(Ri)  < 

<  e  t/e||/e(’)UiO)IU“(H^)  +  1  ~  e  t/e- 


2.  Finite  speed  of  propagation.  We  assume  that  initially,  /e(x,  •,  0)  has  a  compact  support  in 
Rv.  Let  us  first  show  that  /e(x,  •,  t )  remains  compactly  supported.  Indeed,  by  (2.6),  fe(-,  v,  t ) 
and  hence  ue(-,t )  are  uniformly  bounded,  and  therefore  the  contributions  of  Xu(v)(v)  on  the 
RHS  of  (2.2)  are  supported  by  ve[- u^,  where  u =  ||ite(x,  Oil  L°°(RixRi)-  Consequently, 
fe(x,-,t )  given  in  (2.2)  remains  compactly  supported  for  all  t  >  0,  with  support  contained 
in  supp„/e(i,  •,  0)  U[-^oo,  itoo]-  (Note  that  after  an  initial  kinetic  layer  of  order  0(e),  the 
contribution  of  the  initial  data  in  (2.2)  decays  exponentially  fast.  Thereafter,  /e(x,  •,  t)  is  in 
fact  ‘essentially’  supported  in  [— Uqo,  «<»])• 
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With  this  in  mind  we  now  turn  to  prove  the  finite  spatial  speed  of  propagation.  We  shall 
need  a  refined  version  of  estimate  (2.6).  To  this  end  we  first  observe  that  according  to  (2.2) 
-  which  we  rewrite  as 


fe(x,v,t)  =  e  t/ef(x  -  <a(u),u,0)  -!-  (1  -  e  t/e) 


!r=0  e(T  t)/tXue(x-(t-T)a(v),r)(^)dr 

/r=0e(T"t)/^T 


fe(x,v,t)  is  given  by  a  convex  combination  of  fe(x  —  fa(u),u,  0)  and  Xu e(x-(«-T)o(v),r)(w)-  By 
Jensen’s  inequality,  therefore,  we  have  for  any  convex  function,  U(f), 


(2.7)  U(fe(x,v,t))  <  e~£/ef/(/e(x  -  ta(v),v,  0))  +-  f  e(r  t)/'t/(xUl(I-(t-T)a(v)1T)(t'))^. 

Z  Jr— 0 

In  particular,  consider  the  case  U(f)  =  |/|p.  If  we  let  ax  denote  the  maximal  speed  of 
propagation, 


(2.8)  Uoo  =  {max  |a,(u)|  ,  v  e  suppu/e(x, -,  <)}, 

and  B[r]  =  [~r,r\d  C  Ri,  is  the  ball  of  radius  r,  then  (2.7)  implies 

ll/e(-)  Vi  0ll£°°(B[r])  —  e  0)||i,oo(S[r  +  tlloo])  + 

+  -  [  e(T-^/r||x„e(.,r)(v)|Uco(B[r+(t_T)aoo])dr. 

£  Jt-0 

Integrating  the  last  inequality  w.r.t.  v  we  find 

J  WfJ'iVJ)\ !t“,(B[r])^U  —  e  ^  j  !!/«(')  0)||^c*>(fi[r  +  tOoo])<^1,+ 

(2.9) 

+  (1  -  e"'/’)  m3?.  /  ll/.(',».T)lk“(l!|.  +  (l-r)o«,l)<<». 

0<r<t  Ju 

If,  on  the  one  hand,  we  take  the  p-root  of  both  sides  and  let  p  f  oo,  we  obtain 

(2.10)  ||/e(x,u,t)||Loo(B[r]xRv)  <  max{||/e(x,u,0)||L=O(B[r+taoo]xR.),l}, 

in  agreement  with  what  we  had  before,  consult  (2.6).  If  on  the  other  hand,  we  set  p  =  1  in 
(2.9),  we  find  that  the  function  F(t), 

F(t )  =  f  \\M-,v,T)\\Lco(B[T+(t_T)aoo])dv,  0  <  r  <  t, 

J  V 

satisfies 

F(t)  <  e~t/eF( 0)  +  (1  -  e"^‘)  ma jcF(t), 

and  hence  F(t)  <  F( 0).  Thus,  we  have  a  finite  speed  of  propagation  (<  a^)  of  the  uniform 
bound  on  the  moments 

(2.11)  /  ||/e(x,W,*)||i»(B[r])dt7  <  /  ||/e(x,t/,0)||ioo(B[r  +  laG<,])^. 

J  v  Jv 
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In  particular,  the  local  density  is  uniformly  bounded  by  the  initial  data, 

(2.12)  ll“e(*.OIL-(J*-xii+)  <  J\\fc(xtvt0)\\Laa{Fi)dv. 

In  summary,  we  conclude  that  for  initial  data  fe(x,  v ,  0)eL1(Rv;  L°°(RX) )  which  are  compactly 
supported  in  Rv,  the  corresponding  kinetic  solution  fc(x,v,t)  remains  compactly  supported 
in  Rv  and  is  uniformly  bounded  in  Ll(Rv ;  L°°(Rd)),  due  to  a  finite  speed  of  propagation 
<  ««,,  given  by 

Ooo  =  jmax|a,(v)|  |  |u|  <  \\fe(x,  v,  0)||Li(R,.L«(RS))  ,  v  e  suppje(x,  v,  0)J  . 

3.  Monotonicity.  The  signature  function  Xu(v)  is  an  increasing  function  of  u.  Consequently, 
the  fixed-point  iterations  (2.4)  show  that 

(2.13)  fe(x,v,0)  >  ge(x,v,0)  ==>•  fe(x,v,t)  >  ge(x,v,t),  for  all  (i,  u), 

namely,  the  solution  operator  associated  with  the  kinetic  equation  (1.2)  is  monotone.  In 
particular,  if  we  compare  a  given  kinetic  solution  (compactly  supported  in  Rdx  X  R^)  with 
the  steady  state  solutions  Xc<mat.{v)-  if  initially  we  have 

(2-14a)  Xk(v)  <  fc(x,v,0)  or  f,(x,v,0)  <  Xk(v), 

we  obtain 

(2-146)  Xk(v)  <  fc(x,v,t)  or  fe(x,v,t)  <  x#r(w), 

in  agreement  with  (2.6).  And,  since  the  kinetic  solution  operator  is  also  conservative,  the 
Crandall-Tartar  lemma  [5]  implies  the  L1  -contraction  stated  in  (2.5).  In  fact,  at  this  point 
we  can  state  a  little  more,  namely, 

4.  Zd-contraction  revisited.  Taking  into  account  the  finite-speed  of  propagation,  we  can 
repeat  -  along  the  lines  of  Remark  2,  a  localized  version  of  estimate  (2.3)  which  sharpens 
the  L1  -contraction  estimate  (2.5)  into 

(2.15)  II /«(*,«,  0  -  flf«(x,V,OIU>(B[r]xR,)  <  ll/e(*,v,0)  “  3e(^,  0)  ||l>  (B[r  +  tQoo]x  R„)- 

5.  The  various  estimates  quoted  above  indicate  that  after  an  initial  layer  of  order  0(e),  the 
kinetic  solution  asymptotes  to  the  ‘equilibrium  function’,  Xt4x,t)(u)i  where  -  as  will  be  shown 
in  the  next  section,  u(x,  t)  is  the  unique  entropy  solution  of  (1.1). 
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3.  KINETIC  ENTROPY  FUNCTIONS 


Our  analysis  of  the  kinetic  model  (1.2)  hinges  on  the  construction  of  certain  kinetic 
entropy  functions.  A  kinetic  entropy  function  in  this  context  is  a  function,  H(f),  such  that 
as  in  Boltzmann’s  il-Theorem,  any  solution  of  (1.2)  obeys  the  additional  entropy  inequality 

(3.1)  j[dt  +  a(v)  ■  dx)H{fe)dv  <  0. 

We  shall  construct  a  family  of  such  kinetic  entropy  functions  depending  on  extra  fixed 
parameter  k,k  real.  To  this  end,  we  integrate  (1.2a)  against  sgn (/„  —  Xk)  over  the  phase 
space.  Invoking  a  standard  regularization  argument  of  the  signum  function  we  obtain 

(3-2)  f  [dt  +  a(u)  •  dx\\fe  -  Xk\dv  =  --  f  sgn (/,  -  Xk){fc  ~  Xu ,)dv. 

J  v  6  Jv 

Noting  that  the  expression  on  the  right  is  upper-bounded  by 

JvSgn(fe  ~  ~  x»‘)dv  =  ““/[!/«  -  Xk\  +  sgn(/e  -  XfcXxfc  ~  Xu ,)}dv  < 

<  --[/  l/B  -  Xk\dv  -  \ut  -  fc(], 

£  Jv 

we  arrive  at 

THEOREM  3.1.  For  any  solution  f e(.L°° (R* \  L1  {R*  x  Rv ))  of  the  kinetic  model  (1.2),  the 
following  inequality  holds 

(3-3)  J[dt  +  a(v)  ■  dx}\fe  -  Xk\dv  <  -~[J  \  fe  -  Xk\dv  -  | ue  -  fc|]. 

Now,  the  RHS  of  (3.3)  is  clearly  nonpositive  for 

(3-4)  /  |/e  -  Xk\dv  <  --|  /(/e  -  Xk)dv | 

Consequently,  Theorem  3.1  yields 

COROLLARY  3.2.  For  any  k,k  real,  the  following  functions 
(3-5)  H(f.)  =  =  If,  -  Xfcl 

are  kinetic  entropy  functions,  i.e.,  we  have 
(3.6)  j [dt  -f  a(v)  ■  dx]\fe  -  Xk\dv  <  0. 
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Let  us  point  out  that  our  kinetic  entropy  functions,  f/jt(/e),  are  intimately  related  to  the 
entropy  functions  used  by  Krushkov  in  [8].  Indeed,  as  e  |  0  we  expect  (and  later  on  prove) 
that  fe  approaches  \u ■  With  this  in  mind,  the  inequality  (3.6)  turns  into  Krushkov’s  entropy 
inequality  [8] 

d  ci  Q 

(3.7)  —\u  -  k\  +  -5—  [sgn(u  -  k)(Ai(u)  -  A,(fc))]  <  0,  for  any  real  k. 

C/v  i“l 

To  make  this  last  point  more  precise,  we  shall  need  several  lemmata.  We  start  with 

LEMMA  3.3.  Let  feeL00(R( ;  L1(RV\  L°°(Rd)))  be  the  solution  of  the  kinetic  equation  (1.2), 
subject  to  given  initial  data  /e(x,u,  0)  which  are  compactly  supported  in  Rv.  Assume  that 
ue(x,t )  =  fv  fe(x,v,t)dv  satisfies 

(3.8)  ue(x,0)  — >  u0(x)  in  L1(Rd), 
and 


(3.9)  a  subsequence  of  ue(x ,  t)  — >  u(x ,  t)  in 


Then  the  sequence  ue(x,t )  converges  strongly  in  L^c{Rf-,L1{Rdx))  to  u(x,t),  which  is  the 
unique  entropy  solution  of  the  conservation  law  (1-1),  i.e.,  (S. 7)  holds. 


Note:  If  we  take  k  >  ||tt||z,°°(,Rrfx[oiT]))  then  the  entropy  inequality  (3.7)  yields 

Wtu  +  S  ^(u)  -  °- 


i.e.,  ti(x,  t)  is  a  supersolution  of  (1.1);  similarly,  taking  k  <  —  shows  that 

u(x,  t )  is  a  subsolution  of  (1.1).  Hence,  (3.7)  implies  that  u(x,  t )  solves  the  conservation  law 
(1.1). 

To  prove  Lemma  3.3  we  first  prepare 

LEMMA  3.4.  Let  feeL°0(Rf ;  L1(i?v;  L°°(/2^)))  be  the  solution  of  the  kinetic  equation  (1.2), 
subject  to  given  initial  data  fe(x,v, 0)  which  are  compactly  supported  in  Rv.  Then  for  any 
k,  k  real,  we  have 

(3.10)  /  |  fe  -  Xk\dv  -|ue-fc| — >0  in  L^Rf  x  Rdx), 

Jv  e|0 

(3.11)  /  a(v)\fe  -  Xfc| dv  -  sgn(ue  -  k)  f  a(v)(fe  -  Xk)dv — *  0  in  Llloc{Rf  x  Rdx). 

Jv  Jv  ej,0 
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PROOF.  The  vanishing  limit  in  (3.10)  follows  from  the  inequality  (3  3),  for 

(312)  T 

0  <  f  [{[  |/e  -  Xk\dv  -  I uc  -  k\\dxdt  <  -e  [  [  f  [dt  +  a(v)dI\\fe  -  Xki^dxdt — *  0. 

Jo  Jx  Jv  Jo  JxJv  c|0 

To  prove  (3.11)  we  write 

/  I  fe  -  Xk\dv  -  \ue  -  = 

J  V 

(3.13)  =  f  sgn(/e  -  Xk)(fc  ~~  Xk)dv  -  sgn (ue  -  k )  /  (fe  -  Xk)dv  =  /  (/e  -  Xk)s{v)dv. 

J  v  J  v  J  V 

Here,  s(u)  =  s(v;x,  t )  is  the  characteristic  function  given  by 

s(u)  =  sg n(fc(x,v,  t)  -  Xk(v))  -  sgn(ue(x,  t)  -  k). 

Now,  since  s(u)  is  supported  on  the  set 

V  =  {n|sgn(/e  -  Xk)  sgn (uc  -  fc)}, 

and  since 

sgn(/e  -  Xk)  ■  s{v)  =  2,  for  veV, 
we  can  rewrite  (3.13)  in  the  following  form, 

(3.14)  / \fe  -  Xk]dv  -  |t*.  -  fc|  =  /  | fe  -  X/t|sgn(/e  -  Xk)s{v)dv  =  2  f  \fe  -  Xk\dv. 

Jv  JveV  j  vzV 

In  view  of  (3.10),  the  identity  (3.14)  implies 

(3-15)  /  \f,  -  Xk\dv — >0  in  Lj  (R?  x  R*). 

JveV  eiO 

We  conclude  by  noting  that 

/  a(v)\fc  ~  Xk\dv  -  sgn (ue  -  k)  /  a(v)(fe  -  Xk)dv  = 

J  v  Jv 

(3.16)  =  /  a(v)(fc  -  Xk)s{v)dv  =  2  /  a(v)\f*  ~  Xk\dv  <  2ax  f  \fe  -  Xk\dv, 

Jv  JveV  JveV 

and  (3.11)  follows  from  (3.16)  together  with  (3.15).  □ 

Equipped  with  Lemma  3.4  we  turn  to  the 

PROOF  (of  Lemma  3.3).  By  our  assumption  (3.8),  there  is  a  strongly  convergent  subse¬ 
quence  (still  denoted  by)  ue(x,  t)  — >  u(x,  t).  Utilizing  (3.10)  we  obtain 

(3.17)  [ \fe- Xk\dv  =  \ue- k\  =  \u- k\. 

J  V 
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Here  the  overbar  denotes  the  weak  *L°°-limit  of  the  indicated  quantities  after  extraction  of 
appropriate  subsequences,  if  necessary.  (We  note  that  the  existence  of  the  weak  *L°°  limits 
here  and  below  are  justified,  since  in  view  of  (2.11),  fe(x,v,t)  remains  compactly  supported 
in  Rv  and  uniformly  bounded  w.r.t.  e  in  Ll(Rv]  L°°(R*))). 

By  (1.2)  we  have 


and  hence  by  (3.9) 
(3.18) 


ft  ~  X.,  =  -  £[3i  +  a(t>)  •  3,}},  — >  0  in  V, 

c[0 


/  at(v)fcdv  =  /  at(v)Xucdv  =  A;(ue )  =  Afiu). 

J  V  J  V 


This  together  with  (3.11)  gives 


(3.19) 


/  a.(,.')l/.  -  Xk\dv  =  sgn(ue  -k)  a,(u)(/e  -  Xk)dv  = 

J  XJ  J  v 

=  sgn(u  -&)•(/  afiv)fedv  -  /  afiv)xkdv) 

J  X)  J  V 


=  sgn(u  -  &)(A,(ti)  —  At{k)). 

Hence,  in  view  of  (3.17)  and  (3.19),  the  weak  limit  of  (3.6)  recovers  the  entropy  inequality 


(3.20) 


u  -k\  +  Yl  ^~[s8n(w  “  k)(Mu)  -  ^i(*))]  <  0. 


The  above  argument  shows  that  the  strong  limit  of  any  subsequence  of  ue  satisfies  the 
entropy  inequality  (3.20).  Since  the  entropy  solution  of  (1.1)  subject  to  initial  conditions 
(3.8)  is  unique,  we  conclude  that  lim^o  ue(x,  t)  =  u(x,  t)  as  asserted.  □ 

We  now  turn  to  show  that  the  continuum  ‘fluid’  limit  of  the  kinetic  equation  (1.2)  exists 
and  is  governed  by  the  conservation  law  (1.1).  By  Lemma  3.3  it  remains  to  show  that 
ue(x,t)  is  precompact  in  L™c(Rf]  /^(i?*)).  In  this  context  there  is  (by  now)  a  standard 
procedure,  e.g.,  [4],  w7hich  is  based  on  uniform  Bounded  Variation  (BV)  estimate  for  each 
fixed  t,  coupled  with  equicontinuity  (typically,  Lipschitz  continuity),  in  time.  This  brings  us 
to  our  next  lemma  which  states 

LEMMA  3.5.  Assume  that 

\\fc(x,v,0)\\BV{RULl{Rv))=  sup  [ t~  r  [  j  |/e(x  +  Ax,u,0)  -  fc(x,v,0)\dvdx] 

|Ai|^tO  |Ax|  Jx  Jv 

is  bounded  uniformly  in  e.  Then  the  corresponding  kinetic  solution ,  fe(x,v,t),  satisfies 


(3.21) 


-  ll/.(I'U'8)||Bvr(RXxL’(Rv))- 


Moreover,  if  /e(x,  v,  0)eL1(Rv;  BV(R^. ))  are  compactly  supported  in  Rv,  then  we  also  have 
for  t\  ,  t"  '  0, 

(3.22)  || uc(x,t2)  -  Jllii^)  <  1 12  -  fi|  •  ax  ■  ||/«(x>v>0)||bv(**xH(k,))- 


PROOF  S  ince  the  kinetic  model  (1.2)  is  translation  invariant  in  spatial  variables,  we  can 
apply  the  .^-contraction  (2.5)  to  fe(x,v,t)  with  gc  =  fc[x  -f  A x,v,t)  and  obtain  (3.21). 


Integration  of  the  kinetic  equation  (1.2)  over  the  phase  space  yields 

d  ^  d  t 

— ue(x,  t)  +  ^2—j"  at(v)fe(x,  v,  t)dv  =  0, 

and  since  /e(x,  v,  0)eL1(Rv;  BV(R*))  C  Ll(Rv\  L°°(R^.))  is  compactly  supported  in  Rv,  we 
may  use  the  finite  speed  of  propagation  bound  in  Section  2  to  conclude 

f  £2  ^  Q  a 

\\ue(x,t2)  -  ue(x,ii)||Li(fig)  =  /  ||^  —  /  at(u)/t(x)u,r)||Li(R^)dr  < 

Jr=:ti  i_2  OX{  Jv 


<  /  «oo||/e(®,  v,T)\\BV(RixLl(R„))dT. 

J  T~t\^ 

Also,  since  /«.(x,  v,  Q)eL1(Rv:  BV(R*))  C  BV(R *  X  Ll(Rv)),  the  last  inequality  together  with 
(3.21)  imply  the  Lipschitz  continuity  in  time,  (3.22),  which  completes  the  proof.  □ 


Remark.  In  the  course  of  proving  Lemma  3.3,  consult  (3.18),  we  established  only  the  weak 
*L°°  convergence  of  the  spatial  fluxes.  However,  equipped  with  the  BV  setup  of  Lemma  3.5 
we  are  able  to  derive  strong  convergence.  Indeed,  one  may  utilize  the  integral  representation 
(2.2)  to  conclude  that  in  this  case  we  have 

fc(x,v,t)->x u(*,t)(v)  strongly  in  Z,a([0, T]  x  Rdx  x  Rv). 

This  together  with  the  finite  speed  of  propagation  imply 

||  f  ai(v)fe(x,v,t)  -  At{u(x,t)) ||Li([o,T]xRi)  < 

Jv 

<  Ooo  •  || fe{x,V,t)  -  Xu(i,0(V)IUl([0,T]xWx/iv)  * 

elO 

in  contrast  to  the  weak  convergence  stated  in  (3.18).  We  shall  omit  the  details  (consult 
Theorem  3.7  below),  and  we  turn  now  to  summarize  our  results  by  stating  one  following. 


THEOREM  3.6.  Suppose  /e(x,  v,  0)eL1(Rv;  L1  f|  L°°(R*))  such  that 
(3.23)  tte(x,0)  =  J  fe(x,v)0)dv  -»  u0(x)  in  L^R*). 
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Then  the  local  density  of  the  corresponding  kinetic  solution,  ue  =  fv  fe(x,v,t)dv,  converges 
to  the  unique  entropy  solution  of  (1.1),  i.e.,  we  have 

(3.24)  [  f'(x,v,t)dv ->  v(x,t)  in  L°°([0,  T]\  Ll(R*)), 

and  the  entropy  inequality  ( 3.1 )  holds. 

PROOF.  We  begin  by  first  assuming  that  fe(x,v,  0)  is  compactly  supported  in  Ll{Rv ; 
BV(R f)),  uniformly  w.r.t.  e.  By  Theorem  2.1  (consult  (2.12)),  ue(x,t )  are  uniformly 
bounded,  and  by  (3.21)  fhey  have  uniformly  bounded  spatial  variation,  i.e., 

II'U«(I>0I \bV(R£)  ^  Wfe(xyv,t)\\BV(RixLl[Ry))  S  Const. 

Hence  {ue(x,t),0  <  t  <  T}  is  a  bounded  set  in  Llf\BV(Rf)  and  by  Helly’s  theorem  it 
is  therefore  precompact  in  IAJ.K)-  By  (3.22),  ||  ue(x,  t)||r,i(Rg)  is  Lipschitz  continuous  in 
time,  and  by  Cantor  diagonalization  process  of  passing  to  further  subsequence  if  necessary, 

(3.24)  follows.  By  Lemma  3.3  this  completes  the  convergence  proof  for  compactly  supported 

BV  initial  data.  The  general  case  is  justified  by  standard  cutoff  and  BV-regularization  of 
arbitrary  L 1  H^00^*)  initial  data,  consult  [4],  □ 

We  continue  with  a  couple  of  remarks. 

1.  The  kinetic  initial  layer.  We  observe  that  Lemma  3.5  supplies  us  with  an  e-uniform  bound 
on  the  spatial  variation  on  the  microscopic  scale,  (3.21).  The  temporal  variation  (Lipschitz 
continuity),  however,  is  uniformly  bounded  only  on  the  macroscopic  scale,  (3.22).  In  general, 
one  cannot  control  the  temporal  variation  in  the  microscopic  scale  (uniformly  in  e),  unless 
we  can  prevent  the  possibility  of  a  kinetic  initial  layer  in  (1.2).  To  this  end  we  proceed  by 

2.  Preparing  the  kinetic  initial  data.  In  order  to  avoid  a  kinetic  initial  layer,  we  have  to 
bound  §-tfe  uniformly  in  e  and  time,  in  particular  at  t  =  0.  Taking  into  account  the  uni¬ 
form  bound  (in  e  and  t )  of  the  spatial  variation,  (3.21),  it  remains  to  bound  the  nonlinear 
‘interaction’  kernel  on  the  right  of  (1.2),  j(Xu,  —  ft)-  In  particular,  we  therefore  need 

(3.25)  ||  ft(x,v,0)  —  Xu,(*,o)(n)||£,i(Kgxflv)  ^0- 

efO 

Since  by  our  assumption  (3.23)  we  already  have  that 

||XM*.o)W  ~  Xuo(*)(u)|lz,l(fl£xfiv)  =  ~  uo(aOlki(K*)“^00> 


0. 


eJO 


the  requirement  (3.25)  boils  down  to 

(3.26)  ||/e(x,u,0)  —  X«o(*)(t;)||l1(H^x/Jo) 


Thus,  given  the  initial  conditions  u{x,t  =  0)  =  u0(x),  we  have  to  prepare  the  kinetic  initial 
data,  /e(x,u,  0),  such  that  (3.26)  holds.  If  we  prepare  the  kinetic  initial  data  in  such  a 
manner,  then  we  can  derive  explicit  bounds  (uniform  in  time)  on  the  error  between  the 
kinetic  solution  and  the  exact  entropy  solution,  as  told  by 

THEOREM  3.7.  (Error  bound).  We  consider  kinetic  inital  data,  fe(x,  v,  0 )eL1(Rv\  BV{R *)) 
which  are  compactly  supported  in  Rv.  Suppose  we  prepare  the  kinetic  initial  data  so  that 

(3-27)  \\f'{x,v,0)  -  XM*)(v)WlHRZxRv)-^0°- 

Then  the  following  error  bound  holds 

\\fe(x,  v)  0  —  Xtit(x,e)('u)||L1(R| xfiv)  — 

(3.28)  <  2eaoo||/e(x,u,0)|[£,i(fii(;By(Rd))  +  2||/e(x,u,0)  —  Xu0(x)(^)||l1(r^xrv)  *  0- 

e[0 

Consequently,  we  have 

(3.29)  fe(x,  v,  t)  -»  x«(x,t)(u)  strongly  in  L°° (Rf ;  L\Rdx  x  IQ). 


Note.  Preparing  the  kinetic  initial  data  according  to  (3.27)  is  a  strengthened  version  of  our 
assumption  (3.23).  In  this  case,  the  kinetic  distribution  converges  strongly  and  uniformly 
in  time,  to  the  equilibrium  state  Xu,  as  expected.  Also,  all  the  weak  limits  indicated  in 
the  proofs  of  Lemma  3.3  and  3.4  are  in  fact  strong  ones;  in  particular  we  now  have  strong 
convergence  of  the  corresponding  fluxes 

f  at(v)fedv  -»  At(u)  in  L\Rdx  x  Rv)), 

Jv 

compared  with  (3.18). 


PROOF.  Since  the  kinetic  model  (1.2)  is  translation  invariant  in  time,  we  can  apply  the 
Lx -contraction  (2.5)  to  fe(x,v,t)  with  ge  =  fe(x,v,t  -f  At)  and  obtain 

d  d 

(3-30)  ||^/«(*,M)lta(JSXit.)  <  \\g~tfe(x,v,t  =  0)||Li(KgxRv). 

The  kinetic  equation  (1.2a)  enables  us  to  upper  bound  the  RHS  of  (3.30),  namely, 


—fe(x,v,t  =  0)||ii(RgxR.)  < 


(3.31) 


'dt 


1 


<  ||[o(v)  •  dx]fe(x,V,t  =  0)]Ll(RgxRv)  +  ~|| Xu< (x,t=0)  -  fc(x,V,t  =  0)||Ll(RgxRv). 
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The  first  term  on  the  right  of  (3.31)  does  not  exceed 

IIK”) '  dx\fe(x,v,i  =  0)||Li(^xftv)  <  a0o||/e(*,M)||L‘(R.:Bif(flS)); 

the  second  term  is  less  than 

1  2 

-||Xu.(«,*=0)(w)  -  ft(x,v,t  =  OJII^I^XH.)  <  -\\fe(x,  V,  0)  -  X«o  (V)  WlA  (Ag  X  ft.)  • 
Substituting  the  last  two  estimates  into  (3.31)  we  end  up  with 

o 

£\\Qj'fe(x,v>t)\\L'(R*xRv)  < 

(3.32)  <  eaoo||/e(2:,  W,0)||Ll(H^;BU(fl|))  +  2||/e(x,U,0)  —  Xuo(x)(^)||L1(J?fxflv)- 

Finally,  we  use  the  kinetic  equation  (1.2a)  once  more,  obtaining 

\\fe(x,v,t)  —  Xu«(*,t)(v)llz,l(flgxfl*)  ~ 

(3-33)  <  e\\foMx,v,t)\\v(RixRv)  +  c||[a(v)  •  dx]fe(x,v,t) ||li(^xRv)  < 

<  2ea00  Wfe(x,  V,  0)||£,1(h.;BV(/^»  +  2||/e(s,U,0)  -  Xuo(*)(^)|Ui(HgxR„), 
and  (3.28)  follows. 

By  Theorem  3.6  we  also  have  that  ue  —  u  and  consequently  that  x«,  —  Xu  converges 
strongly  and  uniformly  in  time  to  zero,  and  by  adding  this  to  (3.28)  we  obtain  (3.29)  as 
asserted.  □ 

We  note  in  passing  that  the  L1  -contraction  and  the  related  BV  estimates  stated  in  Section 
2  and  Lemma  3.5  are  not  identical  with  the  usual  Ll -contraction  statements  concerning 
viscosity  regularizations  of  entropy  solutions  of  (1.1).  In  fact,  at  any  fixed  time  level,  we 
have 

Wfe  -  5e|U>(**xKv)  >  J  I  J(f,  -  ge)dv\  =  ||u,  -  U/,|| 

By  (3.29),  however,  the  two  statements  coincide  in  the  limit  as  e  [  0, 

Wfe  ~  g,\\ii  (H|xR»)  /  IXu  -  Xw\dvdx  =  ||u  - 

J  xtv 

and  we  recover  the  L1  -contraction  (and  the  corresponding  BV  estimates)  for  entropy  solutions 
of  the  conservation  law  (1.1). 

We  close  this  section  by  calling  attention  to  a  rather  unusual  result  in  the  theory  kinetic 
equations.  Namely,  if  u(x,t )  is  a  smooth  solution  of  the  conservation  law  (1.1),  then  the 
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equilibrium  function  Xu(x,t)(v)  is  an  i1 -solution  of  the  corresponding  kinetic  equation  (1.2). 
That  is 


THEOREM  3.8  (exact  solutions).  Ifu(x,  t)eC  H  ( [0,  T]  x  Rd)  satisfies  the  conservation 
law  (1-1),  then  X*[x,t)(v)  ^  a  kinetic  solution  (1.2)  on  R f  x  [0,T], 

Note.  Theorem  3.8  is  no  longer  valid  when  u(x,t)  contains  shock  discontinuities.  After  the 
formation  of  shock  waves,  the  corresponding  kinetic  solution  has  a  “multivalued”  form,  e.g., 
Xui(x,t)(u)  +  X[u2(x.t),«i3(x,«)]('l’)>  as  in  the  transport  collapse  method  of  Brenier  [2], 


PROOF.  We  have  to  show  that  fe(x,v,t)  =  Xu(x,t)(v )  satisfies  the  kinetic  equation  (1.2a), 
i.e.,  that  for  any  test  function  v,  t ) 

(3.34)  /  /  Xulx,t){v)[dt  +  a(v)  ■  dx\<f>{ x,  v,  t)dxdvdt  -  0. 

Jo  JRi*Rv 

Since  the  integration  in  Rv  is  compactly  supported  (on  [-UocUoo]),  it  is  enough  to  consider 
successively  <f>(x,v,t)  =  <j)(x,  t)  ■  {1,  v,  v2, . . in  which  case  (3.34)  amounts  to  the  equivalent 
conservation  laws 


(3.35) 


d_ 

dt 


up(x, t ) 


d  8 
i=i  t/x* 


/u(x,t) 

vp~1  ai{y)dv 


—  0,  p  =  1,2, . . .  in  V . 


Indeed,  (3.35)  are  the  usual  entropy  equalities  satisfied  by  continuous  solutions  of  (1.1), 
but  violating  (for  p  >  1)  the  Rankine-Hugoniot  conditions  after  the  formation  of  shock 
discontinuities. 


4.  MICROSCOPIC  OSCILLATIONS  AND  COMPENSATED  COMPACTNESS 


In  this  section  we  deal  with  the  one-dimensional  scalar  conservation  law 


(4.1) 


du  dA(u) 
dt  ^  dx 


=  0. 


The  corresponding  underlying  kinetic  equation  reads 


^4'2^  =  ^lXM*.t)(v)  “  /‘KM)],  a(')  =  ^'(O. 

and  we  raise  the  question  of  convergence  of  the  local  ‘particles  density',  ue(x,  t )  =  Jv  fe(x ,  u,  t) 
dv,  towards  the  entropy  solution,  u(x,  t),  of  (4.1).  In  this  section  we  give  an  affirmative  answer 
to  this  question,  which  is  independent  of  compactness  arguments,  i.e.,  the  BV  estimates  used 
in  Lemma  3.4.  Instead,  we  appeal  to  compensated  compactness  arguments,  specifically,  we 
employ  Tartar’s  div-curl  lemma  [13].  In  this  context,  it  is  instructive  to  see  how  oscillations 
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which  persist  on  the  microscopic  scale  are  ‘compensated’  in  a  manner  which  enables  us  to 
pass  to  the  limit  on  the  macroscopic  scale.  We  have 

THEOREM  4.1.  Let  feeL°°(Rf\  L1(RV]  L1  f)  L°°(RX)))  be  the  solution  of  the  kinetic  equa¬ 
tion  (4-2).  Then  ue(x,  t)  =  fv  fc(x ,  v ,  t)dv  converges  strongly  in  LplocfRx  x  Rf ),  p  <  oo,  to  the 
unique  entropy  solution  of  the  nonlinear  conservation  law  (4-1). 

Remark.  The  conservation  law  (4.1)  is  nonlinear  in  the  sense  that  there  exists  no  interval 
on  which  the  flux  A(u )  is  linear,  i.e.,  A'(u )  ^  Const. 

PROOF.  Integration  of  (4.2)  over  the  phase  space  yields 

(4.3a)  dtue  +  dx  f  a(v)fedv  —  0. 

J  V 

The  corresponding  entropy  inequality  reads 

(4.36)  dt  f  | fe  -  Xk\dv  +  dx  [  a(v)\fe  -  Xk\dv  <  0. 

Since  by  (2.11)  the  left-hand  side  of  (4.3b)  lies  in  Murat’s  lemma  [10], [13]  implies 

that  the  negative  measure  on  the  right  of  (4.3b)  lif's  in  a  compact  set  of  Hff{Rx  x  Rf ).  Hence 
we  can  apply  the  div-curl  lemma  [13]  to  the  left-hand  sides  of  (4.3a)  and  (4.3b),  which  gives 

v-c  I  I  fc  -  Xk\dv  -  j  a(v)fcdv  •/]/*-  Xk\dv  = 

(4-4)  _  _  _ 

=  Uc  ■  j  a(u)|/e  -  Xk\dv  -  I  a{y)fcdv  ■  /  \fe  -  Xk\dv. 

J  v  J  v  j  V 

We  recall  that  the  overbar  denotes  the  weak  *Z/°°-limit  of  the  indicated  quantities  after 
extraction  of  appropriate  subsequences,  if  necessary.  Following  [12],  we  can  rewrite  (4.4)  in 
the  equivalent  form 

(4.5)  (ue  -  uc)  •  /  a(v)\fe  -  Xk\dv  =  [ \ft  -  Xk\dv  ■  (  f  a(v)fedv  -  /  a(v)fedv). 

Jv  Jv  Jv  Jv 

Using  (3.10)  and  (3.11),  the  last  equality  is  further  simplified  into 

(4.6)  ( ue  -  uc )  •  sgn(ite  -  k)  f  a(v)(fe  -  Xk)dv  =  \ue  -  k\  ■  (  f  a(v)fedv  -  f  a(v)fedv). 

Jv  Jv  Jv 

We  now  examine  (4.6)  at  an  arbitrary  fixed  location  (x,t);  with  k  =  ue(x,t)  we  find  after 
little  rearrangement 

(4.7)  \ue  -  ue\  ■  (J  a(v)xkdv  -  J  a(v)fedv )  =  0. 
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Of  course,  by  (4.2) 

h  -  Xu .  =  -e[dt  +  a(v)dx\fe  — ♦  0  in  V, 

eiO 

hence 

(4.8)  /  a{y)fcdv  =  [  a(v)xUtdv  =  A(ue ); 

J  V  J  V 

Also,  we  recall  with  k  =  ue(x ,  t )  we  have 

(4.9)  f  a(v)xkdv  =  A{ue). 

J  v 

Inserting  (4.8)  and  (4.9)  into  (4.7)  we  find 

We  -  ue\  ■  (A(ue)  -  A(ue))  =  0. 

This  implies  that 

(4-10)  %)  =  4(ue), 

for  otherwise,  \ue  —  ue\ (x,t)  =  0,  which  in  turn  leads  again  to  (4.10).  Taking  the  weak  limit 
of  (4.2),  we  obtain  with  the  help  of  (4.8)  and  (4.10), 

9  _  d  . 

mu‘ +  aiAM  =  0 

Thus,  (a  subsequence  of)  ue(x ,  t )  converges  to  a  weak  solution  of  the  conservat  r  law  (4.1). 
Moreover,  in  view  of  the  nonlinearity  of  A(u),  equality  (4.10)  implies  that  ue{x,t )  converge 
strongly  in  x  Rf),  1  <  p  <  oo,  consult  Tartar  [13].  Using  this  fact  together  with 

Lemma  3.3  we  conclude  that  ue  converges  strongly  in  Lf^Rx  x  Rt),p  <  oo,  to  the  unique 
entropy  solution  of  (4.1),  as  asserted.  □ 

5.  CONSERVATION  LAWS  WITH  A  SOURCE  TERM 

In  this  section  we  extend  the  above  results  to  inhomogeneous  scalar  conservation  laws 

<3  d  Q 

(5T)  — [u(x,f>]  +  — [Ai(u(x,t))j  =  5(x,f,u),  (x,t)eRdx  x  R+ , 

where  S(x,  t,-)  is  an  L°°(Rd  x  R? ;  C1  )  source  term  satisfying  S(x,  t,  0)  =  0. 

The  corresponding  kinetic  model  equation  reads 

(5.2)  [dt  +  a(v )  •  dx]f'(x,  v,  t)  =  ^[xu.(x,«)(u)  “  Mx>  v>  01  +  S'(x>  l>  V)MX>  v>  0. 
and  is  augmented  with  the  constitutive  relations  (1.2b),  (1.2c). 
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A  unique  kinetic  solution  for  (5.2)  can  be  constructed,  as  before,  by  Banach  fixed  point 
it<  rations  which  yield 

THEOREM  5.1.  The  kinetic  model  (5.2),  (1.2b-c)  is  well-posed  in  L°°(R(\Ll(R^  x  Rv)). 
Moreover ,  if  fe  and  ge  are  two  different  inhomogeneous  kinetic  solutions  of  (5.1),  and  if  we 
let  S'^t)  denote 

(5-3)  =  jrnaxS^x.t,'!;)  |  v  e  supp„/e(x,  v,  t)  (J  supp„5e(x,  v,  t)  J  , 

then  we  have 

(5-4)  max  ||/e(x,u,t)-ge(x,u,t)||Li(fidxfii()  <  e^=o S”(0<1£- \\fe(x,  v,  0) -ge(x,  v, 0)||ll(/^xRv). 

We  shall  only  indicate  the  proof  of  the  L°°(R(;  ff(R^  x  Rv)  stability  stated  in  (5.4).  The 
difference  between  the  kinetic  solutions  fe  —  ge  (with  corresponding  local  densities,  ur(x,  t )  = 
fv  fc(x,v,t)dv  and  we(x,t)  =  fv  ge(x,v,t)dv )  satisfies 

\dt  +  a(v)  ■  dx](fe  -  gt)  =  ^[(x«e(x.t)(u)  -  X».(x,t)(v))  -  (/e  -  ge)\  +  S'(x,t,v)(fe  -  gc). 
Multiplying  this  by  sgn (/„  —  gc )  and  integrating  over  Rv  and  R f  (in  this  order),  we  obtain 

^Il/.(X,U,0  ~  gc{x,V,t)  HlH^xRv)  <  '  H/e(x,  t7,  i)  -  pe(l,V,t)||Ll(^xHv), 

and  (5.4)  follows.  □ 

We  conclude  with  several  remarks  concerning  the  entropy  inequality. 


The  corresponding  inhomogeneous  kinetic  entropy  inequality  now  reads 

(5.5)  [ [dt  +  a(v)  ■  dx}\fe  -  Xk\dv  <  f  S'(x,t,v)\fe  -  Xk\dv. 

Moreover,  by  arguing  along  the  lines  of  the  stability  estimate  (5.4)  we  find  that  for  BV(R*) 
source  terms  we  have 

||/e(x,  v,  T)\\BV(RdxLi(Rv))  <  e^=oS“(t)d‘||/e(x,U,0)||BVqWdxZ,i(fii())+ 

(56)  T  T  , 

+  l/r-^)dT  •  max||5'(x,t,u)l|BK(/?d)dt  •  \\fe(x,v,0)\\Li{RdyR<i). 

This  allows  us  to  keep  the  convergence  statement  of  Theorem  3.6, 

/  /e(x,  v,  t)dv  >u(x,  i),  in  L°°([0,T],  ^(R*)), 

Jv  *10 

in  the  inhomogeneous  case  (5.2).  In  view  of  (5.5),  we  are  also  able  to  recover  the  macroscopic 
‘continuum  limit’  entropy  inequality  for  the  above  limit  u  =  u(x,t),  which  in  this  case 
amounts  to 

yt\ u  -  *|  +  £  ^-[sgn(u  "  k)iMu)  ~  ^.(*0)1  < 

<  sgn(it  —  fc)[5(x,  t,  u )  —  S(x,  t,  fc)],  for  any  real  k. 
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